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A quantitative global skin-friction measurement technique is proposed. An oil-film is doped with a luminescent 
molecule and thereby made to fluoresce in order to resolve oil-film thickness. Particle Image Surface Flow 
Visualization (PISFV) is used to resolve the velocity field of the surface of the oil-film. Skin-friction is then 
calculated at location x as t(x) = pCvi^fx, t) /hfx, t), where t) is the velocity field of the surface of the oil- 
film, h is the thickness of the oil-film, t is time, and fi is the dynamic viscosity of the oil. The data collection and 
data analysis procedures are explained, and preliminary experimental skin-friction results for flow over the 
wing of NASA’s Common Research Model are presented. 
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Nomenclature 

coefficient proportional to quantum efficiency of luminescent molecules 

local chord length 

coefficient of friction 

thickness of oil film 

Intensity of luminescent oil-film 

Intensity of the excitation light incident on the oil-film 

Molarity 

Mean Aerodynamic Chord 
arbitrary integer 
static pressure 

volumetric flow rate per unit depth 

streamwise spacing between fringes (FISF technique) 

time 

velocity field of the surface of the oil-film 

location is space 

skin-friction 

dynamic viscosity of oil-film 
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I. Introduction 


( 1 ) 


S kin-friction, or surface shear-stress, constitutes as much 
as fifty percent of the drag on commercial transport 
aircraft at cruise and can cause loadings on structures 
such as wing surfaces, flaps and landing gear. Also, the 
CFD community relies heavily on skin-friction values when 
setting important parameters such as and y'*' — ^parameters 

which when perturbed only slightly can give rise to very 
different flow field predictions. However, skin friction 
remains a difficult quantity to measure experimentally. Many 
point-based skin-friction measurement techniques have been 
fairly well matured, including oil-film interferometry 
(Naughton and Brown (1996), Naughton and Hind (2013)), 
MEMS (Etebari (2008), Kahng et al. (2002)), and liquid 
crystals (Reda et al. (1998)). Relevant to the present article is 
work done by Tianshu Eiu and his colleagues to develop a 
global luminescent oil-film skin-friction meter which they 
used to experimentally determine skin friction lines for 
several juncture flows, delta wings, and NACA airfoils (Eiu 
et al. (1998, 2008)), their work on optical flow algorithms 
(2008), and the work by Mosharov and his colleagues 
developing Particle Image Surface Flow Visualization 
(PISFV) techniques (Mosharov et al. (2006), Mosharov and 
Radchenko (2012)). The current article aims to build upon the 
existing literature to develop a skin-friction meter which can 
produce quality quantitative skin friction data with 
unprecedented resolution — as many as several hundred 
experimentally determined data points per square inch. 


i(x) = 


h{x,t) 


II. THEORY 

It is shown in Massey (1970) that the volumetric flow rate of 
a thin oil-film of height h subjected only to a pressure gradient 
dP /dx and skin-friction r is given by 


rh^ h? dP 
2/r 3/rdx' 


( 2 ) 


In the absence of a pressure gradient or for a sufficiently thin 
oil-film (h « 1) the second term drops out and 


Q 


Tj^ 

2/r' 


( 3 ) 


Using a control volume approach. Tanner and Blows (1976) 
show that the velocity Vf^ of a particle at the surface of the oil- 
film is then given by 


dQ rh 
dh ii' 


( 4 ) 


Rearranging, we observe that skin-friction may be expressed 
as 


SUMMARY OF TECHNIQUE 

The technique detailed herein is summarized as follows. A 
thin luminescent oil-film is applied to the surface of a model. 
The thickness of the oil-film will later be detected by 
measuring the intensity of the emissions from the luminescent 
molecules in the oil — a brighter region corresponds to thicker 
oil. Fluorescent particles are applied to the surface of the oil- 
film. So that the emissions from the luminescent particles and 
the emissions from the fluorescent oil-film are separable via 
optical bandpass filters, the particles emit at a wavelength that 
is different from that emitted by the fluorescent oil-film. The 
tunnel is operated and images of the oil-film/particle-field are 
taken on regular intervals — approximately six frames per 
minute — alternating between an optical filter for viewing 
emissions from the particles and an optical filter for viewing 
emissions from the oil-film. 

It is emphasized: Two sets of data (images) are taken during 
this experiment. One set of images reveals the particle 
locations on the surface of the oil-film. The second set of 
images captures the luminescent intensity of the fluorescent 
emissions from the oil-film. 

The velocity field of the surface of the oil-film is revealed by 
applying optical flow methodologies to pairs of particle 
images. The luminescent intensity of the emissions from the 
oil-film is directly proportional to the oil-film thickness. 
Therefore, measuring the luminescent intensity of the oil-film 
reveals a global relative thickness distribution of the oil-film. 
Finally, letting /r denote the dynamic viscosity of the oil, h 
denote the thickness of the oil-film, t denote time, and Vf^ 
denote the velocity field of the surface of the oil-film, the 
skin-friction magnitude is calculated at location x and as 


Therefore, in order to determine skin-friction for situations in 
which the pressure gradient is negligible or the oil-film is 
sufficiently thin one needs to know the thickness of the oil- 
film, h, the velocity of the surface of the oil-film, and the 

dynamic viscosity of the oil-film, /r. 

The oil-film is doped with a luminescent molecule so that it 
fluoresces under 450nm light. The luminescent intensity of 
the fluorescence is directly proportional to the film’s 
thickness, h. Eetting Iex(^) denote the intensity of the 
excitation energy being cast upon the oil-film at location x, 
denote the intensity of the light emitted by the 
fluorescent oil-film at location x, and by letting a be a 
constant capturing the quantum efficiency of the dye 
molecule and the concentration of the dye in the oil, 

hm (x, t) = alex ix)h(x, t), (6) 


or 


h(x, t) = 


^em 


(V) 


The velocity of the surface of the oil-film, i;^, is found by 
distributing particles on the surface of the oil-film and 
imaging their motion under testing conditions (wind-on). The 
captured frames are then processed using optical flow 
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algorithms to produce a velocity field for the surface of the 
oil-film. Finally, the dynamic viscosity of the oil is assumed 
to be known, although a temperature-viscosity calibration and 
model-mounted thermocouples would increase accuracy. 


III. Experimental Setup 

The measurement system consists of an excitation source, a 
luminescent oil-film, the fluorescent particles, two optical 
filters, and one camera. A schematic representation of the 
system is shown in Figure 1. 

The viscosity of the silicone oil is chosen to be 10,000cSt. 
This highly viscous oil is chosen because it evolves very 
slowly, allowing the time between successive images to be as 
large as 1 0 to 20 seconds. This large time period provides time 
for the optical filters to alternate between a 550nm bandpass 
filter capturing the particle locations and a 650nm bandpass 
filter capturing the luminescent intensity of the oil-film. The 
highly viscous oil also has the favorable characteristic that it 
remains in the region of interest through wind-tunnel start-up 
and through several different experimental conditions. It must 
be applied as a very thin layer (~100/rm) in order to avoid 
wind-driven ripples at the oil-air interface. This thin film of 
oil also puts the oil-film within the range of thicknesses for 
which the oil-film thickness-luminescent intensity 
relationship is linear. 


CCD Camera 



It is particularly important that the experimenter is able to 
discern between emissions from the particles and emissions 
from the luminescent oil-film. The two substances must 
therefore emit at starkly different wavelengths. Further, there 
is a lower bound on the wavelength of the desired excitation 
source because the experimental community would like to be 
able to pass through windows typical of wind tunnel test 
section the emissions from the excitation source. The 
excitation source was therefore chosen to emit over a 
bandwidth centered at 450nm. The particles were chosen to 
emit at a bandwidth centered around 550nm, and the dye- 
molecule for the oil-film was chosen to emit over a bandwidth 
centered around 650nm. The particles are a PF-33 Palm-Leaf 
Green pigment from Risk Reactor Inc., and the dye for the 
silicone oil is Pyrromethene 650 from Exciton. 


Pyrromethenes are excellent candidates for a dye in silicone 
oil as they are soluble in silicone oil and in toluene, which is 
also a solvent for silicone oil. Optical bandpass interference 
filters are used to separate the emissions from the particles 
from the emissions from the luminescent oil-film (and vice 
versa). These filters are from Melles Griot, are centered at 
550nm and 650nm and have full-width half-max bandwidths 
of 40nm each. The camera is a Photometries CoolSnap K4 
4MP camera, chosen largely due to its availability. The 
excitation source is four 1650mW LED lamps from LED 
Supply. The LEDs emit over a bandwidth that is fairly narrow 
and centered at 450nm. The camera and lights were mounted 
in the plenum above the test section. 

IV. EXPERIMENTAL PROCEEDURE 

This section is designed to outline the procedure one must 
walk through in order to collect skin-friction data using the 
method detailed herein. 

PREPARATION OF OIL, DYE, AND PARTICLES 
Pyrromethene laser dye is added to 10,000cSt silicone oil to 
create a luminescent oil that has a dye molarity of about 
10“^M. This is achieved as follows: 

An empty 50mL jar is placed on a balance and its mass is 
recorded. While on the balance, about three mL of 10,000cSt 
silicone oil is added to the jar. The new mass is recorded. 

The oil must be doped with the Pyrromethene molecules, and 
while Pyrromethene is soluble in silicone oil, it is much more 
readily dissolved into toluene. Approximately 1/lOOOg of 
Pyrromethene is weighed out and dissolved into toluene in a 
25mL volumetric flask. Some of this solution is pipetted into 
a 5mL volumetric flask, and these five milliliters of solution 
are rinsed with toluene into the jar containing approximately 
three mL of silicone oil. Toluene is then added in order to 
reduce the viscosity so that an airbrush is able to dispense the 
solution. Importantly, after the oil-film is applied, the toluene 
will evaporate and only the dye and oil will remain on the 
model. 

The particles are applied by stirring them into toluene and 
spraying them through a second airbrush onto the surface of 
the model. A solution of toluene and particles is prepared: 
Approximately 0.00 Ig particles are added to about lOOmL of 
toluene. This ratio is not particularly important as long as the 
researcher monitors the application and distributes the desired 
number of particles per unit area. The desired number of 
particles per unit area depends on the spatial resolution 
required in measurements. Suitably high particle density is 
preferred for optical flow computation to extract high- 
resolution velocity fields (Liu and Shen (2008)). 

FIDUCIALS 

Fiducials are placed on the model in known locations in order 
to enable camera calibration. Also, because the PISFV — LOF 
method measures sub-pixel displacements of particles, image 
pairs need to align perfectly. Several fiducials are used to 
register image pairs with each other so that only the motion of 
the particles relative to the model is detected. 

REFERENCE IMAGES 

During the experiment several images are taken of the 
luminescent oil-film. These images are used to estimate 
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t). In order to prevent unwanted signal from being 
eaptured by the eamera a 650nm optieal bandpass interferenee 
filter is used for this data eolleetion. However, some photons 
that are not emissions from the lumineseent oil-film make it 
through. One must take eare, then, to subtraet the signal that 
is deteeted by the eamera in the absenee of the lumineseent 
oil-film. Immediately before the lumineseent oil-film is 
applied, a sequenee of images is eaptured and averaged to 
provide an image of the baekground signal that is reeeived by 
the eamera. This image will be subtraeted from subsequent 
images of the lumineseent oil-film taken during the 
experiment. 

The intensity of the exeitation light ean vary spatially in the 
region of interest. This direetly impaets the fluoreseent 
emission from the lumineseent oil-film. Before the 
experiment, a diffuse white surfaee sueh as a sheet of paper is 
plaeed over the region of interest. This surfaee must lie flat 
against the testing surfaee. With the lights arranged just as 
they would be during the wind-tunnel run, a sequenee of 
images is eaptured and averaged. For this, the optieal filters 
are removed, and the integration time and aperture are 
adjusted sueh that the eamera is nowhere saturated. The 
resulting image provides a sealar, by whieh one may 

later seale the lumineseent oil-film images to eorreet for 
spatially varying exeitation energy. 

If in order to ereate different testing eonditions the test artiele 
must be physieally moved during the experiment (sueh as to 
adjust the angle of attaek), then the set of referenee images 
(the oil-film baekground image and the exeitation energy 
distribution image) must be generated for eaeh position. 

APPLICATION OF OIL AND PARTICLES 
In order to avoid eontamination, two separate airbrushes are 
employed — one for the lumineseent oil-film, and one for the 
partieles. The fairly low viseosity oil-Pyrromethene-toluene 
solution is transferred into the eontainer used by the first 
airbrush, and a very light dusting of the solution is applied to 
the testing region. Suffieient time is then allowed to pass so 
that the toluene fully evaporates. The partiele-toluene solution 
is then transferred into the eontainer for the seeond airbrush 
and, making sure that the partieles stay suspended in the 
toluene, the region of interest is dusted with the partieles. The 
researeher should monitor the proeess through an optieal filter 
that permits the green partieles to be readily seen. Time is 
again allowed to pass so that the toluene evaporates. 

DATA COLLECTION 

The wind tunnel is operated. Onee the eonditions in the tunnel 
are on point, data eolleetion may begin. Images are aequired 
at a rate not exeeeding one frame per five seeonds (0.2fps). 
The optieal filter is alternated between 550nm to 650nm for 
eaeh image, produeing two distinet data sets. This is 
illustrated in Figure 2. Even numbered frames are images of 
partieles and odd number frames are images of the 
lumineseent oil-film (or viee versa). After at least 50 images 
have been aequired, the eondition in the wind tunnel may be 
altered and another data set may be aequired. 
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Figure 2. Illustration of image sequence. Optical filters 
alternate so that particle images and luminescent oil-film 
images are acquired in alternating frames. 

V. DATA PROCESSING 

Condueting the experiment generates four data files (sets of 
images): 

1 . Oil-film baekground images 

2. Eight distribution images 

3. Partiele evolution 

4. Oil-film thiekness evolution 

The data sets are used as follows: 

The oil-film baekground image is subtraeted from eaeh of the 
lumineseent oil-film images, leaving behind an image 
whieh is mueh more elosely an image of only the emissions 
from the lumineseent oil-film. 

The light distribution image, denoted /^x, is sealed sueh that 
the peak intensity value is around 1. When one divides a 
lumineseent oil-film image taken during the experiment 
(wind-on) by the light distribution image /^x, the resulting 
image is an image whieh represents the lumineseent intensity 
of a region exeited uniformly by the exeitation souree. This is 
the relative thiekness distribution h = Sinee this is the 

thiekness distribution up to a sealar, the eonstant a is set to 1 
(a = 1). A typieal thiekness distribution image h is shown in 
Figure 3. 
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Relative Oil-Film Thickness Image h 



0 MAC/4 MAC/2 


Figure 3. Relative oil-film thickness map. (e.g. An intensity 
value of 2 corresponds to a region of oil twice as thick as the 
region to which an intensity value of 1 would correspond.) 

Because sub-pixel accuracy is important and because the 
relative positions of the model and the camera tend to change 
during the experiment due to tunnel vibrations, the particle 
images must be registered to each other. The centroids of the 
fiducials are identified and the images are translated into 
alignment. The particle images are then processed through a 
suite of optical flow codes (Liu and Shen, (2008)). The 
resulting velocity field, scaled appropriately, describes the 
motion of the surface of the oil-film in meters per second. A 
typical particle image is shown in Figure 4. 


Typical Particle Image 



Figure 4. Typical Particle Image 


An oil-film relative thickness image h is taken from between 
the two images that were used for the optical flow process. 
That image is registered to the particle images and is used to 
complete the right hand side of Equation 1 . 


VI. Results 

Skin-friction measurements were made on a portion of the 
main wing of a scale model of NASA’s Common Research 
Model (CRM) — a generic transport aircraft which was 
developed for the fourth AIAA Drag Prediction Workshop. 
The model used for the data in this paper is a 3% scale model 
with its wings permanently deflected to match the shape of 
the wings at (full scale) cruise. The vehicle was mounted at 4 
degrees angle of attack, and the tunnel was operated such that 
the Reynolds number based on the mean aerodynamic chord 
was about 675,000. 

The section of the wing used for this study is centered at 75% 
of the half-span from the aircraft’s centerline and has a width 
equal to just over half of the main wing’s mean aerodynamic 
chord. The location of this section is illustrated in Figure 5. 


American Institute of Aeronautics and Astronautics 


5 





Figure 5. Schematic of NASA ’s Common Research Model 
(CRM) detailing where on the model data was taken. 

The model had trip-dots installed about 0.1c from the leading 
edge, where c is the loeal ehord length. The trip dots give way 
to signifieant streamwise vortieity that postpones or prevents 
the onset of separation (shown in preliminary flow- 
visualization data at NASA Ames). The streamwise vortieity 
is readily seen in a developed oil-film downstream of the trip 
dots. This is depieted in Figure 6. All data shown in this artiele 
is taken entirely downstream of the trip dots. 


Close-Up of Oil- Film Downstream of Trip Dots 



3 


Spatial Variation 
in Skin-Friction 


Figure 6. Developed oil-film downstream of trip dots. Skin- 
friction is seen to vary radically in the spanwise direction. 

At the eonditions deseribed in this artiele, the flow was 
everywhere attaehed. Ahead of the trip-dots the spanwise 
eomponent of the flow is signifieant, and aft of the trip-dots 
the flow is largely streamwise, turning outboard only in the 
aft-most five pereent of the ehord. A flow visualization was 
produeed by allowing the flow to shear an oil-film upon 
whieh eonglomerations of fluoreseent pigment were initially 
at rest. The resulting streaklines are visible in Figure 7a. 
Figure 7b is a sehematie of the model and illustrates the 
approximate plaeement of the eamera for the eapture of the 
image in Figure 7a. Skin-frietion lines produeed using optieal 
flow teehniques are shown in Figure 8. 


Skin-frietion data was eolleeted using the Partiele Image 
Surfaee Flow Visualization-Lumineseent Oil-Film (PISFV- 
LOF) teehnique proposed in this paper and by using the 
Fringe-Imaging Skin-Frietion (FISF) teehnique. 
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Skin-Friction Lines From Optical Flow Code 
Background; relative oiFfilm distribution 



Figure 7 a. An image of surface streaklines formed through 
the shearing of fluorescent pigment conglomerations on 
the surface of an oil film. There exists a significant 
spanwise component to the flow just upstream of the trip- 
dots and also along the trailing edge. 



Region Under 
Investigation 


Figure 7b. Diagram showing the perspective from which 
the image in Figure 4. was taken. 



Figure 8. Skin-friction lines produced from a pair of particle 
images using optical flow techniques. An oil-film thickness 
distribution image, h, is in the background. 

The FISF technique is a well-established, minimally intrusive, 
measurement technique and was used to measure skin friction 
on a small section of the wing. The skin friction is found 
through application of a simple expression based on 
lubrication theory relating fringe spacing to skin friction. 
Small drops of oil, of known viscosity, are placed on the wind 
tunnel model surface. The wind tunnel is turned on causing 
the oil to spread, forming a very thin wedge, of nearly linear 
profile. After the flow is turned off, an extended quasi- 
monochromatic light source is used to illuminate the surface. 
Light reflecting from the top surface of the oil will interfere 
with light reflected from the surface of the wind tunnel model 
creating an interference pattern. This pattern is then imaged 
using a digital camera. The distance, AS, between the fringes 
is proportional to the thickness of the oil, which is 
proportional to the skin friction. The skin friction coefficient 
magnitude is calculated by using fringe spacing found from a 
regression algorithm followed by a conversion of the pixel- 
based fringe spacing to physical coordinates using 
photogrammetry. The PC-based application, CXWIN5G, 
developed by Gregory G, Zilliac was used in this study. More 
information about this technique is discussed in the Zilliac 
( 1 996), and past application in Zilliac (20 11) and Ross (20 13). 

Each drop of oil generates one data point. For this 
experiment, only 24 skin friction values were collected during 
one run using the FISF technique. Ideally, more FISF data 
points would be collected for comparison between the two 
techniques discussed. 
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Figure 9a shows a raw image of the sheared droplets of oil. 
The fringes from the destruetive interferenee are evident. 
Figure 9b is a elose-up view of the boxed region from Figure 
9a. The greater values of skin-frietion near the trip dots give 
way to a thinner wedge of oil and therefore a fringe pattern 
that is more spaeed out. The smaller values of skin-frietion 
further downstream give rise to more elosely spaeed fringes. 



Figure 9a. Raw data image taken for fringe-imaging skin- 
friction (FISF) technique. The fringes in the sheared drops of 
oil are visible. The white box is the region cropped for further 
study in Figure 9b. 



Figure 9b. The distance, 
As, between the 
destructive interference 
bands is proportional to 
the thickness of the oil, 
which is proportional to 
the skin-friction. This is 
evident in the image to 
the left, where the 
destructive interference 
bands are far apart in 
the region of high skin- 
friction (A), but closer 
together in the region of 
lower skin-friction (B). 


The data generated by the PISFV-LOF teehnique is a relative 
skin-frietion field. This data is therefore sealed to eome as 
elose as possible (in a least-squares sense) to FISF-measured 
skin-frietion values. Figure 10 depiets the Q-values 
determined by the PISFV-LOF teehnique and the loeations of 
the 24 points at whieh FISF data was taken. 


another mid-ehord loeation. In general the diserepaney seems 
to be around 2% to 15%. 

There exists uneertainty in eaeh teehnique. The viseosity of 
the oil varies by about 1% for eaeh degree Fahrenheit ehange 
in temperature. This translates direetly through Equation 1 to 
a 1% ehange in skin-frietion values per degree Fahrenheit. As 
this partieular tunnel draws ambient air from the room, 
models ean ehange temperature by as mueh as 20 degrees 
Fahrenheit during a single run, and sinee for eaeh set of 
measurements (FISF and PISFV-LOF) temperature was not 
aeeounted for, the diserepaneies shown in Figure 10 are 
within the expeeted range of values. Uneertainties also are 
introdueed beeause of the large spatial gradients in skin- 
frietion magnitude. Near the trip dots the skin-frietion values 
deeay steeply as one moves aft. A slight imperfeetion in the 
measurement of a data-point loeation therefore has a fairly 
signifieant effeet on the aeeuraey of the assoeiated skin- 
frietion value. 


Percent Error of determined by 
PISFV - Luminescent Oil-Film Data 
Relative to 



0 MAC/4 MAC/2 


Figure 10. Map showing Cf -values from the PISFV- 
Luminescent Oil-Film technique together with the percent by 
which those values differ from the Cf values obtained using 
the FISF technique. 


After being sealed to a best fit of the FISF data, the PISFV- 
LOF data shows eneouraging eloseness to the skin-frietion 
values determined at the 24 diserete FISF loeations. The 
pereent by whieh the PISFV-LOF data is different from the 
FISF data is shown in Figure 10. The diserepaney is as high 
as 26% at one loeation near mid-ehord, but as low as 1.6% at 


Streamwise seetions showing the skin-frietion profiles 
measured by eaeh teehnique are presented in Figures 11a and 
11b. Figure 11a shows the regions over whieh the PISFV- 
LOF data was averaged in the spanwise direetion. The FISF 
data points whieh lie in these regions are also plotted with the 
streamwise seetions in Figure 11b. 
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Figure 11a. Regions over which PISFV-LOF data was 
averaged to produce the streamwise sections shown in Figure 
11b. 
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Figure 11b. Streamwise sections PISFV-LOF generated Cf- 
values were averaged in the spanwise direction over the 
regions shown in Figure 11a. The FISF generated Cf -values 
which fell into those regions are also plotted in this figure. 

VII. Conclusion 

The Particle Image Surface Flow Visualization with 
Luminescent Oil-Film (PISFV-LOF) skin-friction 
measurement technique shows promise. The preliminary data 
set shows reasonable trends in skin-friction coefficient 



Region of Spanwise Average in Figure 11b 



magnitudes, and the percent deviations from fringe-imaging 
skin-friction values are fairly low for preliminary data using 
a newly developed technique. Increased accuracy can likely 
be achieved by taking care to account for temperature 
fluctuations of the oil-film. The technique currently produces 
a relative skin-friction field. In order to resolve nominal skin- 
friction values from the technique an accurate oil-film 
thickness-luminescent intensity calibration is necessary. 
Further work could focus on more accurately measuring the 
oil-film thickness through luminescent intensities. 
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